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Abstract 

The result of a search for gaugino pair production with a trilepton signature 
is reinterpreted in the framework of minimal supergravity (mSUGRA) with 
imparity violation via leptonic A Yukawa couplings. The search used 95 pb _1 
of pp collisions at yfs = 1.8 TeV recorded by the D0 detector at the Fermilab 
Tevatron. A large domain of the mSUGRA parameter space is excluded for 
Am, A122 > 10~ 4 . 
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Typeset using REVTpjX 



Supersymmetry (SUSY) is one of the possible extensions of the standard model (SM). For 
each SM particle there is a hypothesized supersymmetric partner with spin differing by 1/2- 
integer. Most searches for supersymmetric particles assume conservation of -R-parity, Rp, a 
multiplicative quantum number defined as (— \} 3Nb+Nl+2S , where Nb is the baryon number, 
Nl is the lepton number, and S is the spin quantum number JT|. However, SUSY does not 
require imparity conservation. In particular, the lightest supersymmetric particle (LSP) can 
decay into a purely leptonic state due to the presence of an R p - and ^-violating term in 
the supersymmetric potential, X^LiLjE^ , where Lj and E\. are isodoublet and isosinglet 
supersymmetric lepton fields, respectively (the superscript C indicates charge conjugation). 
The indices i, j, k run over the three lepton generations and the potential is antisymmetric for 
the indices % and j. Current upper limits on .R-parity violating SUSY Yukawa couplings, A^, 
are of the order of ~ 10 -2 |J. If these couplings are not vanishingly small, an enhancement 
is expected in the number of produced multilepton events. 

In this paper, we reinterpret the result of a previous search by the D0 collaboration 
for gaugino pair production in multilepton channels ||. We use the minimal low-energy 
supergravity model j|,[5| (mSUGRA) as a starting point, and add non-vanishing Ay* cou- 
plings. The mSUGRA model has four continuous parameters and one discrete parameter: 
mo — the universal scalar mass, mi/2 — the universal gaugino mass, Aq — the common 
trilinear interaction term, tan/3 — the ratio of the vacuum expectation values of the two 
Higgs fields, and the sign of /i — the Higgsino mass parameter. The mass spectrum of the 
SUSY partners at the electroweak scale and their decay branching ratios are obtained from 
the above parameters by solving a set of renormalization group equations using the program 
ISAJET 0. Present limits on the Xijk Yukawa couplings [[| imply that this mass spectrum is 
the same as for the case of conserved .R-parity. In this analysis we consider only parameter 
regions with a neutralino (%°) as LSP. 

The CDF and D0 collaborations have previously reported on searches for .R-parity vi- 
olation in the di-electron+jets channels (7]||. They assumed an Rp- and AT/,- violating su- 
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per symmetric potential term \' i j k Q i LjD^ , where Qi and are isodoublet and isosinglet 
super symmetric quark fields, respectively. Some regions of the mSUGRA parameter space 
are excluded by non-observation of SUSY or Higgs particles at the CERN e + e~ collider 
(LEP2): the present limit on the mass of the lightest neutral SUSY Higgs boson (88.3 
GeV P|) implies that tan/3 < 2 is excluded, independent of the other parameters. At higher 



tan/3, part of the parameter space is excluded by the lower limit on the Xi mass [ TOJ obtained 
assuming .R-parity violation through A couplings. 

The event selection and background estimations used in this work are discussed in the 
above-mentioned D0 search ||. Four different final states were considered: eee, ee/i, e/i/i, 
and /i/i/i, requiring at least three electrons, two electrons and a muon, two muons and an 
electron, or three muons, in the respective channels. No acceptable events were found. The 
result is summarized in Table |. The corresponding selection criteria (including the triggers) 
are detailed in Ref. 0. We consider these selection criteria adequate for the present analysis. 

Our search is most sensitive to decays with highest electron and muon multiplicity, i.e., 
those with no r lepton among the decay products of the LSP. The detection efficiency 
is highest, especially for the case of Am, when electrons dominate. On the other hand, 
couplings A133 and A233 correspond to decays with least sensitivity, because the number of 
r leptons is highest. We limit ourselves to the three extreme cases: Am, Am and A233. 

We generate Monte Carlo (MC) events with all possible production and decay modes 
of SUSY particles assuming the mSUGRA model using isajet H with .R-parity viola- 
tion [|ll];|Il] • We apply the same selection criteria as used in || to these generated events, 
and calculate all signal efficiencies. 

Detector response is modeled using a parameterized, fast, particle- level simulation of 
isolated electrons, photons, and both isolated and non-isolated muons. The model contains 
jet reconstruction and a simulation of the missing transverse energy in an event. Lepton 
acceptance criteria include the loss of electrons in the region between the central and end 
cryostats of the calorimeter (1.2 < \r)\ < 1.4), and a lookup table of the muon efficiency 
as a function of r\ and [|T3| , p~4|| , where r\ and 4> are the pseudorapidity and the azimuthal 



angle of the lepton, respectively. The parameters of the program are tuned so that the total 
acceptance, e total ; and the shapes of the missing transverse energy distributions and charged 
lepton rj, 4> and transverse energy distributions agree with detailed simulation based on 



GEANT [p~5| , [T6|1 . The total acceptance includes the geometrical acceptance, efficiency factors 
for the trigger, track reconstruction, and lepton identification. It depends mainly on the 
type of coupling and on the value of m^j. In the vicinity of the exclusion contour, the 
typical values are 20%, 10%, and 0.3% for Am, A122, and A233, respectively. e total decreases 
with decreasing mi/2, mainly because the masses of the gauginos decrease and the energies 
of their decay products fall below the detection threshold. 

Our 95% C.L. exclusion contours are based on a Bayesian approach |p7||18| . For each 



point in the (7%, 7711/2) plane, we calculate a 95% C.L. upper limit on the cross section. The 
excluded region is determined from the intersection of this surface with the corresponding 
cross section predicted by isajet. In this calculation, we use as input the total integrated 
luminosities, and the uncertainties in the numbers of background events (cf. Table |) and 
in e total . The latter includes the statistical error, an overall 10% systematic error in the 
MC simulation, and the error on efficiency factors for the trigger, track reconstruction, and 
lepton identification, determined through independent measurements described in Ref. [^J. 
Their values are between 10% and 20%, and depend on the event category (and therefore 
on the Xijk coupling) and to a lesser extent on event kinematics (e.g., on supersymmetric 
particle masses). Finally, we include a 10% uncertainty on the theoretical cross section, due 
to e.g., the choice of parton distribution function. 

Figures [1] through |] show, respectively, the exclusion regions in the (7%, W4/2) plane 
for the three chosen couplings, for tan/3 = 5 and 10, and for both signs of //. Since the 
characteristics of SUSY signatures at hadron colliders are rather insensitive to values of 
A |Tj|, we have fixed the value of A to zero. The dashed line indicates the limit of 



our sensitivity in mi/2 for the least favorable case, i.e., for the coupling of A233, where 
£ totai < io~ 4 . The exclusion regions correspond to the spaces below the solid lines labelled 
with the coupling types, and above the higher of the dashed line and the dash-dotted curves 
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specifying the numerical values of A. In the regions beyond the dash-dotted curves, the 
average decay length of the LSP calculated for the value of the coupling indicated on the 
curve, is less than 1 cm. Since efficiency studies for high impact parameter tracks have not 
been done, we conservatively restrict the present study to decay lengths less than 1 cm. 
Thus, for example, the region between curves labelled with A121 and 10~ 3 is excluded if 
A121 > 10~ 3 . The shaded areas indicate the regions where there is no electroweak symmetry 
breaking or where the LSP is not the lightest neutralino. Finally, we also show limits 
corresponding to the present lower limit on the Xi mass (dotted line), which exclude the 
regions below. The wiggles on the A233 curves are due to statistical fluctuations and to the 
10 GeV spacing between neighboring m points used to calculate the curves. 

In conclusion, we have reinterpreted the result of a search for trilepton events in terms 
of possible i?-parity violation in decays of the LSP. We have found that a large domain of 
mSUGRA parameter space can be excluded, provided that f?-parity breaking is achieved 
by lepton-number non- conservation with A121 or A122 couplings greater than « 10 -4 . The 
region of sensitivity extends beyond that presently excluded by LEP experiments P,|TU|]. For 



A233, where our experiment is least sensitive, only a very limited domain of parameter space 
can be excluded, and this region is already excluded by LEP. The excluded values of mi/2 
depend mainly on the type of coupling, and much less on the values of other parameters. 
In particular, the excluded region is slightly larger for \i >0 than for fi <0, and is almost 
independent of tan/3. 

We thank the staffs at Fermilab and at collaborating institutions for contributions to 
this work, and acknowledge support from the Department of Energy and National Sci- 
ence Foundation (USA), Commissariat a L'Energie Atomique and CNRS/Institut National 
de Physique Nucleaire et de Physique des Particules (France), Ministry for Science and 
Technology and Ministry for Atomic Energy (Russia), CAPES and CNPq (Brazil), De- 
partments of Atomic Energy and Science and Education (India), Colciencias (Colombia), 
CONACyT (Mexico), Ministry of Education and KOSEF (Korea), CONICET and UBACyT 
(Argentina), A. P. Sloan Foundation, and the Humboldt Foundation. 



TABLES 



Event categories 


eee 


ee/j. 






Ant (P b_1 ) 


98.7 ±5.2 


98.7 ±5.2 


93.1 ±4.9 


78.3 ±4.1 


Observed events 














Background events 


0.34 ± 0.07 


0.61 ±0.36 


0.11 ±0.04 


0.20 ± 0.04 



TABLE I. The result of the search for a trilepton signature at D0 [3]. 
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FIGURES 
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FIG. 1. Exclusion contours at 95% CL. limits for tan/3 = 5, fx <0, for the case of finite A121, 
A122 and A233 couplings. For the explanation of the different curves, see the text. 
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FIG. 2. Exclusion contours at 95% C.L. limits for tan/3 = 5, fx >0, for the case of finite Am, 
A122 and A233 couplings. 
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FIG. 3. Exclusion contours at 95% C.L. limits for tan/3 = 10, fi <0, for the case of finite Am, 
A122 and A233 couplings. 
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FIG. 4. Exclusion contours at 95% C.L. limits for tan/3 = 10, fi >0, for the case of finite Am, 
A122 and A233 couplings. 
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